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Abstract. The theory of holographic algorithms, which are polynomial time algorithms for 
certain combinatorial counting problems, yields insight into the hierarchy of complexity classes. 
In particular, the theory produces algebraic tests for a problem to be in the class P. In this 
article we streamline the implementation of holographic algorithms by eliminating one of the 
steps in the construction procedure, and generalize their applicability to new signatures. Instead 
£SJ of matchgates, which are weighted graph fragments that replace vertices of a natural bipartite 

, | graph Tp associated to a problem P, our approach uses only a natural number-of-edges by 

O ! number-of-edges matrix associated to Tp. An easy-to-compute multiple of its Pfaffian is the 

number of solutions to the counting problem. This simplification improves our understanding of 
the applicability of holographic algorithms, indicates a more geometric approach to complexity 
classes, and facilitates practical implementations. The generalized applicability arises because 
our approach allows for new algebraic tests that are different from the "Grassmann-Pliicker 
identities" used up until now. Natural problems treatable by these new methods have been 
previously considered in a different context, and we present one such example. 



U 
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q 1. Introduction 



In |18 | I19 [ 120 1 121 [ 122. 23J L. Valiant introduced matchgates and holographic algorithms, in order 
to prove the existence of polynomial time algorithms for counting and sum-of-products problems 
that naively appear to have exponential complexity. Such algorithms have been studied in depth 
f^. and further developed by J. Cai et. al. [D EJ EE U El El IE] • 

The algorithms work as follows: suppose the problem P is to count the number of satisfying 
assignments to a collection of Boolean variables x\,... ,x m subject to clauses c\, ... ,Cp. The 
problem P defines a bipartite graph Tp = (V, U, E), with vertex sets V = {x±, . . . , x m } = {xi} 
and U = {ci, . . . , c p } = {c s } and there is an edge (£, s) € E iff xi appears in c s . Holographic 
algorithms apply if the coordinates of the clauses and variables, expressed as tensors, satisfy a 
collection of polynomial equations called matchgate identities, possibly, in fact usually, after a 
change of basis. In the matchgate approach each vertex of Tp is replaced by a weighted graph 
fragment called a matchgate to form a new weighted graph rWp), such that the weighted sum of 
perfect matchings of Tc2(p) equals the number of satisfying assignments to P. If Tq(p) is planar, 
or more generally Pfaffian, the weighted sum of perfect matchings of Tr^p) can be computed 
in time polynomial in \E\ using the FKT algorithm |13} [T7] . FKT defines a sign- altered skew- 
symmetric adjacency matrix X of Ffi(p) whose Pfaffian equals the weighted sum of matchings 
of TQrpy Our approach is related directly to the graph Tp, computing the Pfaffian of a natural 
\E\ x \E\ matrix associated to Tp. We replace FKT with an edge ordering defined by a plane 
curve as described in Section [6} Evaluating the Pfaffian takes polynomial time. 

Equivalently, the number of satisfying assignments to P is the result of the pairing of a vector 
G G C 2 ' B| formed as the tensor product of "local" data representing the variables and a vector R 
in the dual vector space, the tensor product of "local" data concerning the clauses (see Section 



4.1). The Valiant-Cai formulation of holographic algorithms can be summarized as 



[1) ^satisfying assignments of P = (G, R) = weighted sum of perfect matchings of T^p). 



Date: April 2, 2009. 



1 



2 



J.M. LANDSBERG, JASON MORTON AND SERGUEI NORINE 



In this article we give a new construction which eliminates the need to construct matchgates. 
We associate constants a = chg, = (3r (depending only on the number of each type of vertex) 
and \E\ x |£|-skew symmetric matrices z = Zq, y = Ur directly to G, R, without the construction 
of matchgates, to obtain the equality: 

(2) ^satisfying assignments of P = (G, R) = a/3Pfaff(J + y); 

see Examples |4.2.4 and |4.2.5| The constants and matrices are essentially just components of the 



vectors G, R. The algorithm complexity is dominated by evaluating the Pfafnan. 

The key to our approach is that a vector satisfies the matchgate identities iff it is a vector of 
sub-Pfaffians of some skew-symmetric matrix, and that the pairing of two such vectors can be 
reduced to calculating a Pfaffian of a new matrix constructed from the original two. A similar 
phenomenon holds in great generality discussed in Appendix ^7J A simple example is the set of 
vectors of sub-minors of an arbitrary rectangular matrix. We describe an example of such an 
implementation in Section |5j 

The starting point of our investigations was the observation that the matchgate identities 
come from classical geometric objects called spinors. The results in this article do not require 
any reference to spinors to either state or prove, and for the convenience of the reader not familiar 
with them we have eliminated all mention of them except for this paragraph and an Appendix 
(§|7]), included for the interested reader. However further results, such as our characterization of 
1-realizable signatures [T3], do require use of the representation theory of the spin groups. 

2. Counting problems as tensor contractions 

For brevity we continue to restrict to problems P counting the number of satisfying assign- 
ments of Boolean variables X{ subject to clauses c s (such as #Pl-Mon-NAE-SAT in Example 



3.1). Following e.g., [2] express P in terms of a tensor contraction diagrammed by a planar bi- 
partite graph Tp = (V, U, E) as above (see Figure [T]), together with the data of tensors Gi = G Xi 
and R s = R Cs attached at each vertex Xj £ V and c s € U. Gi will record that X{ is or 1 and 
R s will record that the clause c s is satisfied. Let n = \E\ be the number of edges in Tp. 

For each edge e = (i, s) S E define a 2-dimensional vector space A e with basis a e | , a e |i- Say 
Xi has degree di and is joined to Cj 15 . . . , Cj d . Let Ei denote the set of edges incident to X{ and 
associate to each Xi the tensor 

(3) Gi :=a^ Sh | <g> • • ■ <g> a ijSjd j + a^-ji® • ■ • <g> a*,^ |i e A4 := A iSji ® • • ■ ® A iSj ^ 

= <S> a e i + ® o e h 

The tensor Gi represents that either xi is true (all l's) or false (all O's). It is called a generator 
and in the matchgates literature is denoted by the vector (1, 0, . . . , 0, 1) corresponding to a 
lexicographic basis of Aj = ® e ^EiA e . This vector is called its signature. We use notation 
emphasizing the tensor product structure of the vector space A\ = C 2dl , and will use the word 
signature to refer to the tensor expression of Gi. 

Next define a tensor associated to each clause c s representing that c s is satisfied. Let A* e be 
the dual space to A e with dual basis a e | ,a e |i. Let E s denote the set of edges incident to c s . 
For example, if c s has degree d s and is "not all equal" (NAE), then the corresponding tensor 
(called a recognizer) associated to it is 

( 4 ) R s-= Yl a i,si\ei®---®<Xi,s da \e ds = Yl % a e\e e 

(ei,...,e ds )^(O r ..,0),(l,...,l) (e 1 ,...,e ds )^(0,...,0),(l,...,l) ^ * 

Now consider G:=(&iGi and R:=(& S R S respectively elements of the vector spaces A:=(& e A e 
and A* := (g> e ^4*. Then the number of satisfying assignments to P is (G,R) where (•,•) : 
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A x A* — > C is the pairing of dual vector spaces. At this point we have merely exchanged our 
original counting problem for the computation of a pairing in vector spaces of dimension 2^1 



A 1 

Variable Vertices V 
Edges E 



Clause Vertices U 




As 

^q^j G = G 1 (g>G 2 0G 3 e A 1 (g>A 2 



R = R!®R 2 e Al®A%®Al®Al 



Figure 1. A bipartite graph T diagrams a tensor contraction, representing an 
exponential sum of products such as counting the satisfying assignments of a 
satisfiability problem. Boxes denote clauses, circles denote variables. Each clause 
or variable corresponds to a tensor lying in the indicated vector space; e.g. R\ G 
A*<g)j4.|g)A3. Instead of replacing each vertex with a matchgate, our construction 
defines an n x n matrix, where n is the number of edges in the problem graph. 
The Pfaffian of this matrix, times a constant depending on the number of each 
type of variable and clause, is the number of satisfying assignments. 



3. Local conditions and change of basis 

In order to be able to construct the matchgates corresponding to the Xj, c a , there are local 
conditions that need to be satisfied; the algebraic equations placed on the Gi, R s are called the 
Grassmann-Pliicker identities (or Matchgate Identities in this context). See, e.g., Theorem 7.2 of 
|3] for an explicit expression of the equations, which are originally due to Chevalley in the 1950's 
[9]. These identities ensure that a tensor T representing a variable or clause can be written as 
a vector of sub-Pfaffians of some matrix. From the matchgates point of view, these equations 
are necessary and sufficient conditions for the existence of graph fragments that can replace the 
vertices of Tp to form a new weighted graph Tq(p) such that the weighted perfect matching 
polynomial of ^n(p) equals (G, R). 

Expressed in the basis most natural for a problem, a clause or variable tensor may fail to 
satisfy the Grassmann-Pliicker identities. However it may do so under a change of basis; e.g. in 



Example 3.1 we replace the basis (True, False) with (True+False, False— True). Such a change 
of basis will not change the value of the pairing Ax A* — > C as long as we make the corresponding 
dual change of basis in the dual vector space — but of course this may cause the tensors in the 
dual space to fail to satisfy the identities. Thus one needs a change of basis that works for both 
generators and recognizers. In this article, as in almost all existing applications of the theory, 
we only consider changes of bases in the individual AgS, and we will perform the exact same 
change of basis in each such, although neither restriction is a priori necessary for the theory. 

3.1. Example. In #Mon-3-NAE-SAT, we are given a Boolean formula in conjunctive normal 
form where each clause has exactly three literals, and all are either positive or negative (no mixed 
negations). A clause is satisfied if it contains at least one true and one false literal. The counting 
problem asks how many satisfying truth assignments to the variables exist. The generator tensor 
Gi corresponding to a variable vertex Xj is (§. The recognizer tensor corresponding to a NAE 
clause R s is Q and in our case we will have d s = 3 for all s. 
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Let To be the basis change, the same in each A e , sending a e o *— > a e |o + a e|i an d a e |i a e|o — a e|i 
which induces the basis change a e | i— > ^(a e |o + a e|i) an d a e h •— ► 5(a e |o — a e|i) m ^t- This basis 
is denoted b2 in |23j . Applying To, we obtain 

T o(a it s n \o® ■ ■ ■ ® ai,s ld .\o+ai,s n \i®---® a>i,s id ,|i) = 2 ^ °i,sn|ei® " " " ® a i,s id . |e d . 

{(6i,...,e d .)|Ee^=0(mod2)} 

In the matchgates literature this tensor is denoted by the vector (2, 0, 2, 0, . . . , 2, 0, 2) (assuming 
the number of incident edges is even). We also have 

V(ei,62,£3)^(0A0),(l,l,l) 

= 6a Ml |o®ai iS2 |o®a M3 |o - 2(a Ml | ®a M2 | 1 (g>a iiS3 |i + a i ^ 1 \i®a i>S2 \o®a. i ^\i + a i!Sl \ 1 ®OL i ^\i<g l a^ S3 \Q) 
or denoted by its coefficients, (6, 0, 0, 0, —2, —2, —2, 0). 

4. Holographic algorithms without matchgates 

Though we do not use matchgates, in our approach the matchgate identities still must be 
satisfied under a change of basis as above. Our purpose is to make G, R expressible as vectors 
of sub-Pfaffians of some skew-symmetric matrices. While this appears to be a global condition, 
it can be accomplished locally. Let u, v be s x s and t x t matrices, and form a block diagonal 
(s + 1) x (s + 1) matrix from them. The vector of sub-Pfaffians of the new block diagonal matrix 
in C 2 " ' can be obtained by taking the 2 s x 2* matrix corresponding to the product of the 
(column) vector of sub-Pfaffians of u with the (row) vector of sub-Pfaffians of v , and writing the 
matrix as a vector in C 2S+ \ The analogous statement holds for block diagonal matrices built 
out of an arbitrary number of smaller matrices. Thus if each Gi,R s is a vector of Pfaffians, 



the corresponding G, R will be so too; see Proposition 4.1.2. Theorem 4.2.2| below shows how 
realizing G, R as vectors of sub-Pfaffians aids one in computing the pairing (G, R) indicated by 
<§. 

In all this there is a problem of signs that we have not yet discussed. The problem arises 
because if we order the Xi and c s , there are two natural types of orders of the vector spaces 
in the tensor products of the Ai s , one grouping i's and one grouping s's. The block- diagonal 
discussion above cannot be simultaneously applied to both orderings at once. We explain this 
problem in detail and how to overcome it in ^6] 

4.1. The complement pairing and representing G and R as vectors of sub-Pfaffians. 

Assume we have a problem expressed as above and have constructed tensors G, R such that in 
some change of basis their component tensors Gi,R s satisfy the Grassmann-Plucker identities. 
For the purposes of exposition, we will assume the total number of edges is even. See the 
discussion in the Appendix ^7] for the case of an odd number of edges. 

To compute {G,R), we will represent G and R as vectors of sub-Pfaffians. For an n x n 
skew-symmetric matrix z, the vector of sub-Pfaffians sPf (z) lies in a vector space of dimension 
C 2 , where the coordinates are labeled by subsets I C [n], and 

(sPf(z)) 7 = Pfaff(z/) 

where zi is the submatrix of z including only the rows and columns in the set /. Letting 
I c = [n] \ I, similarly define sPf G C 2 " by 

(sPf v (z))/ = Pfaff(z jc ). 
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The vector spaces A, A* come equipped with un-ordered bases induced from the bases of the 
A e . These bases do not have a canonical identification with subsets of (1, . . . , n) but do have a 
convenient choice of identification after making a choice of edge ordering. After an ordering E 
of the edges has been chosen we obtain ordered bases of A, A*. To obtain the convenient choice 
of identification for A, identify the vector corresponding to I = (ij.,... , i2p) with the element 
with l's in the ii, . . . , %i v slots and zeros elsewhere, so, e.g. 1 = corresponds to (0, . . . , 0), 
/ = (1, . . . , 2n) corresponds to (1, . . . , 1). Reverse the correspondence for A*. 

For later use, we remark that with these identifications, as long as the first (resp. last) entry 
of G (resp. R) is non-zero, we may rescale to normalize them to be one. (If say, e.g., the first 
entry of G is zero but the last is not, and last entry of R is non-zero, we can just reverse the 
identifications and proceed.) Note that the first and last choices of entries are independent of 
the edge ordering, but if necessary, to get the first and last entries non-zero, we simply take a 
less convenient choice of identification. (See ^7] for an explanation of this freedom.) As long as 
this is done consistently it will not produce any problems. 

In the rest of this section we assume that the local problem has been solved, i.e., that Grassmann- 
Pliicker identities hold for all the Gi , R s possibly after a change of basis. We also assume for 
brevity that the Gi and R s are symmetric, i.e. that Gi = sPf(xj) = sPf(7r(xj)) for any permuta- 
tion 7r on the edges incident on Gf, this covers most problems of interest. For the more general 
case when the variables or clauses are not symmetric, and we need to be more careful about 
defining Eq and Er, see Section [fj] and the Appendix fj8] 

Definition 4.1.1. Call an edge order such that edges incident on each Xi £ V (resp. c s £ U) 
are adjacent a generator order (resp. recognizer order) and denote such by Eq (resp. Er). 

Proposition 4.1.2. Suppose P is a counting problem as above, Eq and Er are respectively 
generator and recognizer orders. If for all Xi £ V there exists Zi £ Mat^xdi such that sPf (zj) = 
Gi under the Eq identihcation, and similarly, there exists y s £ Matd s xd s for c s and R s with 
sPf v (y s ) = R s , then there exists z,y £ Mat\ E \ x \ E \ such that 

sPf (z) = G under the Eq identihcation and 

sPf (y) = R under the Er identihcation; 

z, y are just given by stacking the component matrices Zi, y s block-diagonally. 

As the proposition suggests, a difficulty appears when we try to find an order E that works 
for both generators and recognizers. 

Definition 4.1.3. An order E is valid if there exists skew-symmetric matrices z, y such that 
sPf(z) = G and sPf v (y) = R under the E identification. 

Thus if an order is valid 

(G,i?) = ^sPf / (^)sPf/ c (y) 

and in the next subsection we will see how to evaluate the right hand side in polynomial time. 
Then in ^6] we prove that if Tp is planar, there is always a valid ordering. 

4.2. Evaluating the complementary pairing of vectors of sub-Pfaffians. Let n be even. 
For an even set I C [n], define cr(I) = and define sgn(7) = (— l) f7 ( / )+l / l/ 2 . Proofs of the 

following lemma can be found in |1 6|, p. 110] and |11[ p. 141]. 
Lemma 4.2.1. Let z and y be skew-symmetric n x n matrices. Then 

n 

Pfaff(z + y) = Yl Yl sgn(/)Pfaff/(z)Pfaff / c(y) 

p=0 /C[n],|/|=2p 
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To use Lemma 4.2.1 to compute inner products we need to adjust one of the matrices to 

Let z be an 



correct the signs. For a matrix z define a matrix z by setting z\ 



n x n skew-symmetric matrix. Then for every even I C [n], 

Pfaffj(S) = sgn(J)Pfaff j(z). 
This is because for odd \I\, both sides are zero. For |/| = 2p, p = 1, . . . , |_§J, 

Pfaff/(z) = (-l) i i+ i 2+ 1 • • • (-l) i2 p- 1+i2 p +1 Pf a ff / (z) = sgn(/)Pfaff 7 (z). 
Thus we have the following Theorem. 

Theorem 4.2.2. Let z,y be skew-symmetric n x n matrices. Then 

(sPf(z),sPf v (y)) =Pfaff(z + y). 

In Section [6] we show that if Tp is planar there is an easily computable valid edge ordering 
E. Our result may be summarized as follows: 

Theorem 4.2.3. Let P be a problem admitting a matchgate formulation T = (V, U, E) (e.g., a 
satisfiability problem as in the first paragraph) such that 

(1) There exists a change of basis in C 2 such that all the Gi,R s satisfy the Grassmann- 
Pliicker identities (i.e., all Gi and R s are simultaneously realizeable) with complementary 
indexing. 

(2) There exists a valid edge order (e.g. ifT is any planar bipartite graph) 

Normalize vr(G) (resp. t(R)) so that the first (resp. last) entry is one, say we need to divide 
by a,(3 respectively (i.e. a = n« a i where Gi = ajsPf(xj) and similarly for (3). Consider 
skew symmetric matrices x,y where is the entry of (the normalized) tt(G) corresponding to 
I = (hj) an d Vj is the entry of (the normalized) t(R) corresponding to I c = (i,j). Then the 
number of satisfying assignments to P is given by a/?Pfaff(5 + y). 

Example 4.2.4. Figure [2] shows an example of #Pl-Mon-3-NAE-SAT, with an edge order given 
by a path through the graph. The corresponding matrix, z + y is below. In a generator order, 





^7 

(«0 




Figure 2. Example Egg and the term S = (1, 2)(3, 6)(4, 5)(7, 8)(9, 12)(10, 11) 



in the Pfafhan (which has no crossings). 

each variable corresponds to a f° ) block. In a recognizer order, each clause corresponds to 
a 3 x 3 block with —1/3 above the diagonal. Sign flips z i— ► z occur in a checkerboard pattern 
with the diagonal nipped; here no flips occur. We pick up a factor of ^ for each clause and 2 
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for each variable, so a = 2 6 , (3 = (|j) 4 , and a/3Pfaff(5 + y) = 26 satisfying assignments. 
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Example 4.2.5. Another #Pl-Mon-3-NAE-SAT example which is not read-twice and its z + y 
matrix are shown in Figure [3] The central variable has a submatrix which is again ones above 
the diagonal and also contributes 2 to a, so a = 2 5 , [3 = (t^) 4 - Four sign changes are necessary 
in z. The result is a/3Pfaff (z + y) = 14 satisfying assignments. 
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Figure 3. Another #Pl-Mon-3-NAE-SAT example and its z + y matrix. 



5. Beyond Pfaffians 

As mentioned in the introduction, the key to our approach is that the pairing of a vector in a 
vector space of dimension 2™ with a vector in its dual space can be accomplished by evaluating 
a Pfaffian if both vectors are vectors of Pfaffians of some skew-symmetric matrix. This type 
of simplification occurs in many other situations as explained in the Appendix §[7] below. One 
simple such is that if the vector space is of dimension (^) and the vectors that are to be paired 
are vectors of minors of some k x (n — k) matrix. Then the pairing can be done by computing 
the determinant of an easily constructed auxiliary (n — k x n — k) or {k x fc)-matrix, so if k is 
on the order of l^j there is a spectacular savings. Explicitly, for k x £ matrices z and y, with 
G = sDet(z) and R = sDet(y), 

(G,R) = det(Id+z T y). 
Here is an example that exploits this situation. 

Example 5.0.6. Given a graph G and an arbitrary orientation of E(G), the incidence matrix 
B = {bl) ve v{G),e&E{G) is a \ V{G)\ x \E(G)\ matrix defined by 

1 if is the initial vertex of e, 
— 1 if v is the terminal vertex of e, 
otherwise. 

For W C V(G) and F C E(G), with \W\ = \F\, let A.w,f(B) denote the corresponding minor 
of B. Let sDet(£?) = (1, e B, . . . , A\y t p(B), . . .) denote the vector of minors of B. 
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A rooted spanning forest of G is a pair {H,W), where W C V(G), H is a spanning acyclic 
subgraph of G, and every component of H contains exactly one vertex of W. The minor A^f(B) 
equals to ±1 if (G\f, V(G) — W) is a rooted spanning forest, and Aw,f(B) = 0, otherwise. 
(See [10] for a proof of a generalization of this statement to weighted graphs.) Therefore, the 
value of the pairing 

(sDet(B*),sDet(B))= £ ]T (A WtF (B)) 2 

WCV(G) FCV(G) 

is equal to the number of rooted spanning forests of G. It is shown [10] that this value can be 
computed efficiently by the Cauchy-Binet formula: 

]T (A W ,F( J B)) 2 = det(Id+ J B* J B), 

WCV(G) FCV(G) 

where Id is a x |.E(G)| identity matrix. 

From our point of view the result outlined in this example is an instance of the above fact 
that the pairing of vectors in the Grassmannian and its dual can be computed efficiently. (The 
Grassmannian can be locally parametrized by vectors of minors of matrices.) The above efficient 
algorithm for enumerating rooted spanning forests is surprising in the same sense as many 
holographic algorithms are: A closely related problem of enumerating spanning forests of a 
graph is #P-hard [12] . 

6. Edge ordering and sign 

Throughout this section we assume the local problem has been solved and we only need a 
valid edge order. We do not require symmetric signatures. 

Given an order E we would like to know if it is valid. Say Eq, Ep are generator and recognizer 
orders so that there exist skew-symmetric matrices z, y such that with respect to these orders 
G = sPf(z), R = sPf(y). Let tt,t G ©i^i respectively be the permutations such that tt(Eg) = E 
and t(Er) = E. Then for all J C [n], Pfaff,^ J){^(^)) = sgn(7r| j)Pfaff j{z) and similarly for r, so 
up to signs we have what we want. Valid orderings yield it, t which preserve sub-Pfaffian signs. 

We now describe one type of valid ordering for planar graphs, called a C -ordering. For any 
planar bipartite graph Tp, a plane curve C intersecting every edge once corresponds to a non- 
self-intersecting Eulerian cycle in the dual of Tp and can be computed in 0( | .£7| ) time. Fix 
such a C, an orientation and a starting point for C, and let E c be the order in which the 
resulting path crosses the edges of Tp. Define Eq to be the generator order chosen so that the 
permutation 7r : Eq — > E G is lexicographically minimal. In particular, Eq agrees with Eg on 
the edges incident to any fixed generator in V. For example, the generator order on Figure [2] is 
1,2,3,6,4,5,7,8,9,12,10,11. Define £g similarly 

To show that E c is valid we will need another characterization of the sub-Pfaffians and the 
notion of crossing number. Let S = {(ei, e[), . . . , (e^, e' k )} be a partition of an ordered set I, 
with |/| = 2k, into unordered pairs. Assume, for convenience, that e r < e' r for 1 < r < k. Define 
the crossing number cr(S) of S as 

cr(S) = #{(r, s) | e r < e s < e' r < e s }. 

Note that cr(S') can be interpreted geometrically as follows. If the elements of I are arranged on 
a circle in order and the pairs of elements corresponding to pairs in S are joined by straight-line 
edges, then cr(5) is the number of crossings in the resulting geometric graph (see Figure [2^b) ) . 
When the order E on / is unclear from context we write cr(S, E), instead of cr(S'). 

For / C E(T), denote by T/ the subgraph of T induced by /. Let y(Tj) be the set of 
pairings S = {(ei, e[), . . . , (e^, e' k )} of I such that edges in each pair share a vertex in the set 
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V of generators. In other words, (e^e'j) £ S implies there exists j £ V, s,t £ U such that 
e i = (j; s )i e i = (i)^)- I n what follows we focus on generators, the corresponding statements for 
recognizers will be clear. 

Proposition 6.0.7. Let T be a bipartite graph and let Eq be a generator edge order. Assume 



the hypotheses of Proposition 4.1.2 are satisfied with z the skew-symmetric \E\ x \E\ matrix 



such that sPf(z) = G with the order Eq. Let I C [n] = E. Then 

Gj = Pfaff,(*) = (-ir (S) *S 

where z s is the product il^^es z e it ^- 
Proof. Let cr(S) denote the permutation 

<r(S) = (ei e[ e 2 e' 2 ... e k e' k ). 

By |15[ p. 91] or direct verification, sgn(<r(5)) = (— l) 01 ^. Therefore, for a skew-symmetric 
matrix z one has 

Pfeff,(*) = ^(-1)"^, 

where zs '■= z eie ^ . . . z eke i and the sum is taken over the set 5? of all partitions of / into pairs. 

We need to show that the terms zs, S £ r y \ J^(Tj) are zero. Note that for a nonzero term, 
there must be an even number of edges in the restriction to each variable. If S contains a pair 
with split ends (xjC s , x^Ct), i^k, then z$ = 0. □ 



6.0.7 



The analogous statement to Proposition 6.0.7 holds for recognizers. We can now prove the 
following Lemma. 

Lemma 6.0.8. Let P be a problem as above such that all the associated Gi,R s satisfy the 
Grassmann-Pliicker relations under some change of basis, Tp is planar and let E c be a C- 
ordering. If tt, z are defined as above, then sPf(7r(z)) = 7r(sPf(z)). 

Proof. It suffices to show that for any / C E(T) and any partition S £ S^iTi) of /, the signs of 
the term corresponding to S in Pfaff/(,z) and Pfaff 7r (/)(-7r(z)) are identical. By Proposition 
this is equivalent to showing that 

( 5 ) (_l)cr(W = JJ ( _ 1) cr(S[.,Eg) j 

where the left hand side of (|5| is the sign of the term corresponding to S appearing in Pfaff 7r m(7r(^)), 
and the right hand side is the sign of the term corresponding to S in Pfaff i(z), as 

Pfaffj(z) = n pfaff /|>)- 

Here S\ x and I\ x denote the restriction to the edges incident to x of S and I, respectively. 

A stronger equality, namely cr(S,E c ) = X^eV cr (^U) Eq), holds. The curve C determining 
E separates V from U. To exploit the geometric intuition presented above, we replace each 
vertex in x £ V by a small circle and join the ends of edges in / on this circle by line seg- 
ments corresponding to pairs in S\ x . The total number of crossings in the resulting graph is 
YlxcV cr (^'lx) Eq) = YlxeV CI (S\xi E C ): as and E c coincide on the set of edges incident to a 
fixed x £ V. On the other hand, a pair {r, s} is counted in cr(S", , E ), if and only if the curves 
with ends on C corresponding to e r U e' r and e s U e' s cross. 

In other words, we are considering restrictions of (the union of e r and e' r ) and (the union e s 
and e' s ) to the region of the plane bounded by C containing V. □ 
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It follows from Lemma 6.0.8 and a symmetric statement for r that E c is valid. 



Example 6.0.9. An example is given in Figure |4j There, the curves composed of edges (3 
and 5), and (4 and 6) cross, and that shows that the permutation (3 5 4 6) is odd. The edges 
corresponding to, say, (3 and 5) and (2 and 7) don't cross, and the permutation (3 5 2 7) is 
even. In the example, 

5 = {{1,9},{2,7},{3,5},{4,6},{8,10}}. 

The term corresponding to S in Pfaff 7r m(7r(x)) is (— l^xi^a^jXs^a^gXs.iOi as ct(S,E c ) = 2. 
The term in Pfaff (x) is a product of —£3^x4,6 an d — xi ) 9X2, 7^8, 10) which are the terms in Pfaffians 
of blocks corresponding to x and x' , respectively. 
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7. Appendix: Spinors and holographic algorithms 

The Grassmann-Pliicker identities are the defining equations for the spinor varieties (set of 
pure spinors). These equations date back at least to Chevalley in the 1950's [S]- The spinor 
varieties, of which there are two (isomorphic to each other) for each n, §+,§_, respectively live 
in A even C n =: 5+, and A ocW C n =: 5_. The parity condition corresponds to requiring that G, R 
both be either in 5+ or 5_. If n is odd then <S+, 5_ are dual vector spaces to one another, and if 
n is even, each is self-dual. It is this self-duality that leads to the simplification of the exposition 
with n is even - the discussion for n odd is given below. 
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They admit a cover by Zariski open subsets where each subset in e.g. §>+ is covered by a map 
of the form 

(6) : A 2 C n -► ®jA 2j C n = A even C n 

(7) x h+ (Pfaffj(x)) 

as I C (1, ... , 2n) runs over the subsets of even cardinality (and by convention Pfaff0(x) = 1). 

The identification 5+ ~ ^even^2n - g no ^. canon i ca j. \y e can get different identifications by 
composing <p with the action of the Weyl group. The Weyl group action assures that some "less 
convenient" map will have first entry nonzero for G, R as mentioned in §4.1| 

The map ^ is a special case of a natural map to the "big cell" in a rational homogeneous vari- 
ety and the potential generalizations to holographic algorithms mentioned to in the introduction 
would correspond to replacing S + by a Lagrangian Grassmannian or an ordinary Grassmannian 
of /c-planes in a n-dimensional space. More generally, if V is a generalized G(n)-cominuscule 
module, where n denotes the rank of the semi-simple group G, then the pairing V x V* — » C, 
when restricted to the cone over the closed orbits in V, V* can be computed with 0(n 4 ) arith- 
metic operations, even though the dimension of V is generally exponential in n. 

Much of the exposition could be rephrased more concisely using the language of representation 
theory. For example, the fact that if each Gi lies in a small spinor variety then G = ®Gi lies 
in a spinor variety as well, is a consequence that the tensor product of highest weight vectors 
subgroups with compatible Weyl chambers will be a highest weight vector for the larger group. 
Similarly the map z \—* z has a natural interpretation in terms of an involution on the Clifford 
module structure that 5+ comes equipped with. 

On the other hand C 2 may be viewed as (C 2 )® n and as such, inherits an SZ^C-action. The 
SL2 (C) action corresponds to our change of basis, and what we are trying to do is determine 
which pairs of points can by simultaneously be moved into the spinor varieties in (C 2 )® n and the 
dual space (C 2 *)® n . The convenient basis referred to in the text corresponds to an identification 
that embeds the torus of SL2 diagonally into the torus of Spiri2n so weight vectors map to 
weight vectors. 

To continue the group perspective in complexity theory more generally, one can also view 
the ability to compute the determinant quickly via Gaussian elimination as the consequence of 
the robustness of the action of the group preserving the determinant: whereas above there is a 
subvariety of a huge space (the spinor variety) on which the pairing can be computed quickly, 
and a group SL2 that preserves the pairing - a holographic algorithm can be exploited if the 
pair (G,R) can be moved into the subvariety S + x S + under the action of SL2. In Gaussian 
elimination, for the corresponding subvariety one takes, e.g., the set of upper-triangular matrices, 
and the group preserving the determinant acts on the space of matrices sufficiently robustly that 
any matrix can be moved into this subvariety (and in polynomial time). Contrast this with the 
permanent which is also easy to evaluate on upper-triangular matrices, but the group preserving 
the permanent is not sufficiently robust to send an arbitrary matrix to an upper-triangular 
one. This difference in robustness of group actions might explain the difference between the 
determinant and permanent, as well as why only solutions to certain SAT problems can (so far) 
be counted quickly. 

8. Appendix: Non-symmetric signatures 

Most of the natural examples of holographic algorithms, and, in particular, the examples 
given in this paper, correspond to generator and recognizer signatures Gi and R s which are 
symmetric, that is invariant under permutations of edges incident to the corresponding vertex. 
The assumption that the signatures are symmetric is also convenient for our arguments. If 
the signatures are symmetric, then the generator tensor G can be represented as a vector of 
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sub-Pfaffians in some generator order if and only if it can be represented as such a vector in 
every generator order, and the same holds for recognizer orders. This does not hold for general, 
non-symmetric signatures. We now explain how to deal with non-symmetric signatures. 

It is shown in Section |6] that given a planar curve C, an edge order E c and a generator order 
Eq, the tensor G can be represented as a vector of sub-Pfaffians in E c if and only if it can 
be represented as one in Eq. A similar statement holds for E and the recognizer order E^. 
The edges incident to a given generator are ordered in a clockwise cyclic order in Eq. It is easy 
to verify that only a cyclic, not linear, ordering enters Grassmann-Plucker identities. Thus for 
non-symmetric signatures the following statement holds. 

Theorem 8.0.10. Let P be a problem admitting a matchgate formulation Tp = (V,U,E) 
with r p planar. Let the edges incident to every vertex ofTp be ordered in a clockwise order. 
Assume that there exists a change of basis such that all the Gi, R s satisfy the Grassmann-Plucker 
identities with complementary indexing. Then there exists a valid order, and the number of 
satisfying assignments of P can be found in polynomial time. 

Note that the assumption that the edges (or "wires" ) are ordered in a way that agrees with a 
planar embedding is also used in the matchgate formulation, as the matchgates must be inserted 
in such a way that the resulting graph remains planar. 



